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I. INTRODUCTION 



Superdeformation and clustering in ""^Ca studied with Antisymmetrized Molecular 

Dynamics 

Y. Kanada-En'yo and M. Kimura 
Yukawa Institute for Theoretical Physics, Kyoto University, 
Kyoto 606-8502, Japan 

Deformed states in ^"Ca are investigated with a method of antisymmetrized molecular dynamics. 
Above the spherical ground state, rotational bands arise from a normal deformation and a superdefor- 
mation as well as an oblate deformation. The calculated energy spectra and E2 transition strengths 
in the superdeformed band reasonably agree to the experimental data of the superdeformed band 
starting from the Ojj" state at 5.213 MeV. By the analysis of single-particle orbits, it is found that 
the superdeformed state has particle-hole nature of an 8p-8h configuration. One of new findings is 
parity asymmetric structure with ^■^C-|-^*Si-like clustering in the superdeformed band. We predict 
I that ^^C-f^^Si molecular bands may be built above the superdeformed band due to the excitation 

■ of inter-cluster motion. They are considered to be higher nodal states of the superdeformed state. 
' We also suggest negative-parity bands caused by the parity asymmetric deformation. 
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' In the spectra of ^"Ca, the existence of low-lying deformed states has long been known since 1960s [1]. The 

• coexistence of the spherical and various deformed states in the doubly magic nuclei ^°Ca has attracted a great 
^ [ interest. In the positive parity states below 7 MeV, besides the spherical ground state, the existence of two deformed 

• rotational bands built on the state at 3.35 MeV and the 0^ state at 5.21 MeV is known experimentally. The 

■ structure of these bands has been often discussed in relation to many-particle many-hole states of nuclei in this mass 
■T^j" : region for a long time. Grace and Green [2,3] firstly suggested that the former band is dominantly a ip-Ah state, while 

the latter band with a large deformation is understood as an 8p-8h state. 

The member states of the first K'^=0^ rotational band built on the state are known to be strongly populated 
by a-transfer reactions and to have the large a spectroscopic factors [4-7]. The a-l-'^^Ar cluster structure of this band 
has been studied theoretically based on semi- microscopic cluster model calculations [8-10]. 

The states in the second band have been observed in the experimental work searching for the 8p-8h states 

' with ■^^S(^^C,a)'*°Ca reactions, which were performed by Middleton et al. in 1972 [11]. Due to the strong population 
I in the multi-nucleon transfer data and the strong E2 transitions [1], the 0"'"(5.21 MeV), 2"'"(5.63 MeV) and 4+ (6.54 
Q MeV) have been thought to belong to the superdeformed band with a dominant 8p-8h configuration. Recently, by 
^ using the GAMMASPHEAR array detectors, the level structure of the deformed bands in ^°Ca has been explored 
Ph ' and many excited states up to very high spins have been discovered [12,13]. Consequently, the existence of the 
J> [ superdeformed band of ''°Ca built on the 0j^(5.21 MeV) has been experimentally well established. 
• ■ contrast to the experimental situation, theoretical interpretation of the superdeformed band is not well established 

' and is still complicated. With the symmetry restricted Skyrme-IIartree-Fock(HF) calculations, Zheng et al. studied 
: the energy systematics of nparticle-nhole(np-n/i) states in ^°Ca [14,15]. They suggested an 8p-8h deformed states 
■ - - ' as a local minimum of the energy as a function of n. They also predicted highly deformed 12p-12h states in a high 
excitation energy region. However, the recent mean-field calculations suggest that the 8p-8h configuration becomes 
unstable, especially when the model space is extended. Based on the three-dimensional(3d) HF calculations without 
any assumption of the spatial symmetry, Inakura et al. studied the energy surface of '*'^Ca as the function of the 
quadrupole deformation. They have shown that the local minimum of the superdeformed 8p-8h configuration on the 
energy curve becomes quite shallow or even disappear due to the 7 and octupole deformations depending on the 
effective interaction, though the 12p-12h state and an oblate state remain to be local minima in any of used effective 
interactions. They have pointed out that the extreme softness of the superdeformed 8p-8h configuration against the 
octupole deformation (Y30 and Y31). The instability of the 8p-8h configuration and the stability of the 12p-12h and 
oblate solutions in '^''Ca had been also reported by Leander and Larsson [17] based on the macroscopic-microscopic 
model calculations. Further calculations beyond mean-field theory were performed by Bender et al. with GCM 
calculations in HF+BCS approach [18]. In their results, the superdeformed band is not dominated by a certain np-nh 
configuration but it contains various configurations such as Ap-Ah, 6p-6h, and 8p-8h states due to the strong mixture 
between these configurations. Moreover, shell model calculations were performed for description of the superdeformed 
bands [19]. Although the level structure and quadrupole moments of the superdeformed band are described by the 
spherical shell model with the fixed 8p-8h configurations, however, when the model space is extended to include other 
particle-hole configurations, the quadrupole moments are much underestimated due to the configuration mixing. 
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These facts may imply that a stable solution for the intrinsic state of the superdeformation is hardly obtained in the 
modern calculations with mean-field approaches. It seems to be somehow inconsistent with the experimental energy 
spectra of the superdeformation in which no large deviation from the rigid rotor model is seen [12,13]. Therefore, it 
is natural to expect that there might be some mechanism beyond mean-field approaches to stabilize the superdefor- 
mation. The cluster aspect is one of the essential features for nuclear deformation even in sd-pf shell nuclei together 
with the mean-field aspect. Indeed, it has been already shown that two kinds of nature, i.e. cluster and mean-field 
coexist in the deformed states in sd-shell nuclei [20-22]. We expect that the clustering can play an important role also 
in the formation of the superdeformed state of *°Ca. In case of ^°Ca, we consider that the softness against octupole 
deformations in the superdeformed state suggested by Inakura et al. can be associated with asymmetric di-cluster 
structure, though, to our knowledge, there has been no studies to investigate its cluster aspect. 

We also note that some srf-shell nuclei have the remarkable cluster aspect, molecular resonances, whose properties 
have been investigated by cluster models. Recently, there are theoretical efforts to link the molecular resonances with 
low-lying deformed states with full microscopic calculations. It is suggested that the molecular resonances and the 
low-lying deformed states in srf-shcU arc regarded to be a series of molecular bands. For example, in '^^S system, 
it has been suggested that ^^O-^^O molecular resonances are built on the superdeformed band of ^^S [22] owing to 
the excitation of the inter-cluster motion. Therefore, we consider that molecular resonances can be good probes to 
understand the cluster aspect of the low-lying deformed states. In ''"Ca system, there are experimental implications 
of ^^Si-f-^^C molecular resonances in the low-energy fusion cross section [23] and also elastic scattering data [24], 
where distinct structures have been found in the energy region £'cni=20-30 MeV above the ^^Si-|-^^C threshold. It is 
interesting to look into possible appearance of molecular resonances and its relation to the superdeformation in ^"^Ca. 

In this paper, we study deformed states in ''°Ca with antisymmetrized molecular dynamics(AMD). The method 
of AMD has been proved to be useful in description of shape coexistence as well as cluster structure in srf-shell 
nuclei [20-22,25 29]. One of advantages of this method is that both cluster aspect and mean-field aspect can be 
described within the AMD framework. In order to study the coexistence of deformations in ^°Ca and to investigate 
the corresponding rotational bands, we apply a constraint AMD method. Parity projection, which is essentially 
important to describe the parity asymmetric shape of nuclei, is performed before energy variation, while total-angular- 
momentum projection is done after the energy variation. We also perform superposition of the wave functions with 
different configurations obtained by the constraint AMD to obtain better wave functions and to describe possible 
configuration mixing. Level scheme and E2 transition strengths are calculated from these superposed wave functions. 
In the present work, we adopt the effective nuclear interaction which contains finite-range 2-body and 3-body forces 
proposed in Ref. [30]. This interaction is suitable for describing the nuclear structure properties such as binding 
energies and radii over a wide mass number region from ''He to ''"Ca. Moreover, it is also useful to describe nucleus- 
nucleus potential. These features can not be described by zero-range forces such as the Skyrme forces, and are superior 
points of the effective forces with finite-range three-body terms. We believe that they are essential to represent cluster 
aspect and to investigate properties of many-particle-many-hole states with large deformations. 

This paper is organized as follows. In section II, the formulation of AMD is briefly explained. In Sec. Ill, we show 
the calculated results of the energies, quadrupole transitions and moments, and also the intrinsic structure. We give 
discussions of molecular aspect and single-particle orbits in Sec. IV. Finally, a summary is given in Sec.V. 

II. FORMULATION 

We here briefly explain the formulation of the present calculations. For the details of the AMD formulation for 
nuclear structure studies, the readers are referred to Refs. [20,31,32]. 

A. Intrinsic wave function 

We start from the intrinsic wave function for a A-nucleon system which is expressed by a Slater determinant of 
single-particle wave packets, 

|$AMd) = -^A{<fi,<f2, fA}, (1) 

iPi{rj) = M^j)XiVi, (2) 

where A is the antisymmctrizcr, and (pi is the ith single-particle wave packet consisting of the spatial part the 
spin part Xi and the isospin part r]i. The spatial part of the ith single-particle wave packet is written by a Gaussian 
wave packet located at Xj. 



2 



Mrj) = exp{-u{rj-Xif}, (3) 

Xi = (^+^i)XT + (^-Ci)Xi, (4) 
r]i = proton or neutron. (5) 

Here, the center Xj of Gaussian is a complex number and treated as a independent free parameter for each nucleon. 
The width parameter v takes a common vahic for ah nuclcons, and is chosen to be an optimum value for each nucleus. 
The orientation of the intrinsic spin is expressed by a variational complex parameter ^i, and the isospin function is 
fixed to be up(proton) or down(ncutron). Thus, an AMD wave function is specified by a set of variation parameters, 
Z = {Xi, X2, • • • , X4, ^1, ^2, • ■ ■ ,^a}, which are optimized by the energy variation. As the variational wave function, 
we adopt the parity projected AMD wave function, 

I^^md) = ^^I^amd). (6) 

Thus the parity projection is performed before the variation. This is essentially important to describe the parity 
asymmetric intrinsic structure such as the asymmetric cluster structure as discussed in Refs. [20,33]. 



B. Hamiltonian and constraint 



The Hamiltonian consists of the kinetic energy, nuclear and Coulomb forces. In the present work, we use the 
effective nuclear force consisting of two-body and three-body forces. Then the Hamiltonian is written as, 

i i<j i<j<k 

where ii is a kinetic operator, ■O^^^ and v^^^ are two-body and three-body forces, respectively. Tcm is the kinetic term 
for the center-of-mass motion. In the AMD wave function, the energy of center-of-mass motion is subtracted exactly. 

We apply a constraint AMD method, where we find the minimum energy solution in the AMD model space under 
certain constraints. Since the particle-hole structure of the low-lying state is one of the main interests in the present 
work, we adopt the constraint of a total number of harmonic oscillator(H.O.) quanta, 



^os^E 



4h'u ' 2 



(8) 



where the width parameter v is taken to be the same value as that of the single-particle wave packets in Eq.(3). This 
is the principal quantum number of the spherical harmonic oscillator with the oscillation number co = 2'hv/m{m is 
the nucleon mass). In order to obtain the minimum energy state with the constraint {Nos) = Nqs, where A^os is a 
given number, we perform energy variation for the expectation value £ of the Hamiltonian, 



with respect to the variational parameters Z by the method of frictional cooling [20] under the constraint. We note 
the AMD wave function obtained with the energy variation for positive or negative parity state under the constraint 
{Nos) = as $AMD(iV+) or <^amd{N~), respectively. 



C. superposition 

In the obtained wave function $amd(A^os^''), there exist local minimum solutions in the energy curve as a function 
of Nos- As shown later, those are usually the local minima also with respect to the deformation parameter (3 and 
considered to be approximate intrinsic wave functions for the corresponding deformed bands. In the calculation of 
level scheme and physical obscrvablcs, we choose several intrinsic wave functions <I>AMD(fc) and superpose the spin- 
parity projected wave functions, -pMif ^amd(^)' order to satisfy the orthogonality among the calculated states 
and to obtain better wave functions. Here, k indicates Not^ and stands for the label for the obtained intrinsic wave 
function($AMD (A; = Not^)). The wave function of the nth J"^ state is written as, 
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(10) 



k,K 



where coefBcients are determined by digonalizing the Hamiltonian and norm matrices: 



{PiiK'^iMD{k')\H\Pi,K"^iuD{k")) and 

i^MK' ^AMD i^') I ^MK" ^AMD )) ■ 



(11) 

(12) 




(13) 



If we use an enough number of base wave functions, this procedure of the superposition is equivalent to the generator 
coordinate method with respect to the generator coordinate iVos. However, the number of the basis in the superposition 
is 10 at most because of huge computational time in the numerical calculations of the total angular momentum 
projection in the present work. It means that J2k stands for the summation of discrete basis, and hence, the present 
calculations do not correspond to the complete GCM calculations. 



We use an effective nuclear force which consists of central force, spin-orbit and Coulomb forces. As the central 
force, we use F3B force [30] , which contains a finite-range three-body term. The parameter set (2) of the F3B force is 
adopted. As for the spin-orbit force, we use the same interaction as that adopted in Ref. [30], the two-range Gaussian 
spin-orbit term in G3RS force [34] . In the present calculations, the strength parameter of the spin-orbit force is chosen 
to be {i)uis = ui = —U2 = 2500 MeV or {ii)uis = 1800 MeV. Coulomb force is approximated by the sum of seven 
Gaussians. 

In Fig. 1, we show the binding energies oi Z = N nuclei from *He to *°Ca obtained by using the present interaction 
parameters. The calculations are the simple AMD calculations without constraints. The binding energies calculated 
by variation after parity projection are shown in Fig. 1(a), and the energies with total-angular-momentum projection 
after the variation with parity projection are plotted in Fig. 1(b). We also show the results of Ref. [30] for the 
interaction(ii), where AMD calculations with constraint and superposition were performed. The energies for such 
deformed nuclei as ^Be, ^^C, ^''Ne and ^^Mg are gained by the total-angular-momentum projection. As shown in 
Fig. 1(b), the binding energies of medium sd shell nuclei are overestimated by the interaction(i) due to the stronger 
spin-orbit force than that of the case(ii). 



We apply the constraint AMD method to "'"Ca. The width parameter v for the single-particle wave packets is chosen 
to be as = 0.14 fm~^ which optiniizes the energy of the ground state. We impose the constraint of the total number 
of harmonic oscillator quanta as (iVos) = A^os on the parity projected AMD wave function, and obtain energy curves 
as a function of A^os- With a given constraint value A'os, we obtain a few states with different shapes as local minima. 
A prolate shape and an oblate shape are obtained by randomly choosing initial states in the variational procedure. 
In order to find possible cluster states, we also start the variation from initial states with various di- nuclear structure 
and decrease the A'os- Then we find that a a+^^Ar-like clustering appears as a local minimum solution. 

The energy curves as a function of the oscillator quanta(A'os) are shown in Figs. 2 and 3. In principle, the minimum 
value of A^os for the positive-parity state is 60, which is given by the srf-shell-closed configuration, and it is 61 for 
negative-parity states. The results shown in the figure are calculated with the total-angular-momentum projected 
wave functions -Pmx^amd(-^o^^)' where AT = is chosen to obtain the lowest energy. The energies are also plotted 
as a function of quadrupole deformation /32. Here the deformation parameter /32 is defined by using the sharp edge 
liquid drop relation between Q2 and /J2 as follows; 



III. RESULTS 



A. Effective nuclear force 



B. Shape coexistence 
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FIG. 1. The binding energies oi Z = N nuclei from ''He to '"'Ca obtained by the simple AMD calculations without constraints. 
The binding energies calculated by variation after parity projection but no total-angular momentum projection are shown in the 
upper panel(a), and the energies obtained by total-angular-momentum projection after the variation with parity projection are 
plotted in the lower panel(b). Squares are the experimental binding energies. Filled circles and crosses indicate the results by 
using the interaction (i) and (ii), respectively. The results for the interaction(ii) in Ref. [30] of AMD calculations with constraint 
and superposition are shown by triangles with dotted line. 
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where the expectation vahies are calcuhited for the intrinsic states $ ^j^.j^ ( A^os ) , and z is chosen to be an approximately 
symmetric axis. Here we use the fixed radius parameter, i?, which is related to the nuclear mass A = 40 according to 
the formula R = IflA^I^ fm. 

We first mention the results of positive-parity states with the interaction(i) (Figs. 2(a) and (b)). We find three 
states with different shapes as the local minima on the energy surface. The lowest state in the = 0^ curve is 
an almost spherical state at iVos = 62. There are two local minimum solutions in the energy c;urvcs, which almost 
degenerate to each other, at 7~8 MeV higher energy than the lowest solution. One is an oblate state at A'^os ~ 65 and 
the other is a large prolate deformation at N^^ ~ 70 with /32 ~ 0.52. As shown later, the former contains a dominant 
Ap-Ah configuration, and the latter is dominated by a '8p-8/i' configuration and corresponds to the superdcformed 
state. The appearance of the oblate state is consistent with the macroscopic- microscopic calculation [17] and also 
the 3d-HF calculations [16]. On the other hand, the calculated energy curve for the prolate states differs from the 
macroscopic-microscopic calculation by Lcander and Larsson [17] where there is no local minimum for the 8p-8h 
configuration but is a highly deformed 12p-12h state. The a-f-^^Ar-like cluster state appears in A'os > 69 region, 
but it is much higher than the superdcformed state in the case of interaction(i). With a smaller N^g value of the 
constraint, the a-cluster is absorbed into the '"'^Ar core and the a+'^^'Ar-like state changes into the normal prolate 
state during the energy variation. In the 2+ states, the energy curve for the prolate states is more gentle than the 0"*" 
curve at the small deformation region, < (32 < 0.2. 

Next we look into the results obtained with the interaction(ii) (Figs. 2(c) and (d)), which has the weaker spin-orbit 
force than the interaction(i). Although there exist the oblate state and the superdeformed state as local minima 
as same as in the results(i), the excitation energies of these states are much higher than those in the case of the 
interaction(i). The excitation energy of the superdeformed state is about 20 MeV, which largely overestimates the 
experimental energy, 5.21 MeV, of the bandhead of the superdeformed band. On the other hand, the energy of 
a-clustcr state does not so much change from (i) to (ii) because the effect of the spin-orbit force is smaller than other 
excited states. As a result, a-cluster state almost degenerates with the superdeformed state and does not vanishes 
even in the smaller A^os region(as A'os > 66) in the results(ii) than the case(i). Comparing Figs. 2(c) and (d), we find 
that the total number of H.O. quanta has a good correspondence with the quadrupole deformation j32 in the TVqs > 64 
region for the prolate deformation. On the other hand, in the small Aqs region, all the solutions for A'os ^ 63 have 
almost spherical shapes with /?2 ~ 0. It means that different configurations with the same 02 value are obtained as 
optimum solutions for the given A'os. 

We illustrate the density distribution of the intrinsic state ^AMY){Not^) in Fig. 4 and Fig. 5. A remarkable point is 
that parity asymmetric structure appears in the prolate states. Especially, one of striking features is the prominent 
cluster structure like ^^C+^^Si in the large prolate deformation. Comparing the results of (i) and (ii), we remark that 
the intrinsic structure for a given A^os is almost the same as each other. It means that the structure of coexisting 
shapes is not sensitive to the choice of either interaction (i) or (ii) except for the a-cluster state. 

In contrast, the excitation energies strongly depend on the strength of the spin-orbit force as shown in the comparison 
between the results with a stronger spin-orbit force (i) and with a weaker one(ii). Particularly, the oblate state and 
the superdeformed state appear at the low excitation energies in the results with the interaction(i), but they exist at 
much higher excitation energies with the interaction(ii). Considering the Ap-Ah configuration in the oblate shape and 
the ^^C-^^Si clustering in the superdeformed state, it is easily understood that the excitation energies of these states 
decrease with the stronger spin-orbit force(i) than with the weaker one(ii). Although the systematics of the binding 
energy in the sd-shell region seems to be reproduced better by the force(ii) than the force(i), however, the excitation 
enc;rgy of the superdeformed state is reasonably reproduced with the forc;c(i) while they are largely overestimated by 
the force(ii). The force(i) also reproduces the excitation energy of the lowest negative-parity state(J'^ = 3~) rather 
well than the force(ii) as shown later. Therefore, we examine mainly the results(i), because our interest is in excited 
states of ^"Ca. especially in the properties of the superdeformed state. We should note again that the structure of the 
superdeformed state is stable in both cases (i) and (ii), and most of the intra-band properties of the superdeformed 
band does not so much depend on the choice of the interaction(i) or (ii) except for the relative energy to the ground 
state and other rotational bands. 

Now we turn to negative-parity states. In Fig. 3, we show the energies of = 1~ (-Pm^$amd(-^os ')) ^'^'^ 

i^MK^AMui^o^ ''))' obtained with the interaction(i), as functions of (a) A'os and (b) l32- It is found that the 
constraint on (Nqs) acts as the constraint on deformation /?2 in negative-parity states as well as in the positive-parity 
ones. In the small deformation region, the lowest state is the = 3~ which corresponds to a Ip-lh state. The 
difference between the minimum energy of the 3~ state and that of the 0+ state is 4.7 MeV, which well corresponds to 
the experimental excitation energy 3.74 MeV of the lowest 3~ state. We here comment that the interaction (ii) gives 
a much larger energy difference as 6.9 MeV than the interaction (i) and fails to reproduce the experimental excitation 
energy of the 3~ state. As A^qs becomes large, the J'^ = 1~ state becomes lower than the J'^ = 3~ state and a 
K"^ = 0~ rotational band with a prolate intrinsic state is built on this state. We illustrate the density distributions of 
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FIG. 2. The energy curves for the positive-parity states obtained by the constrained AMD calculations. The constraint of 
the total number of harmonic oscillator quanta, (iVos) = A^os, is imposed on the parity projected AMD wave function. The 



energies of the spin-parity projected states, 



MO ^AMdI^'o' 



arc plotted (a)(c) as a function of A'os and (b)(d) as a function 



of deformation /32 of the states: ^amd(-^o^')- The left panels((a) and (b)) show the results with the interaction (i), and the 
right panels((c) and (d)) are for the results with the interaction(ii). 



the intrinsic states for the negative parity states, <I>amd(^os ) in Fig.4(g)-(i). For each Nos, it is found that intrinsic 

structure of i>AMD(-Afos~^) for the negative-parity state is quite similar to that of the prolate solution $amd (-^vis^' ) 
for the positive-parity state. In particular, the negative-parity state with Nos > 68 has the parity-asymmetric shape 
because of the clustering which is regarded as the -'^^0+^^ Si- like structure as well as the case of positive parity state. 
As a result, the largely deformed state forms the K"^ = 0~ rotational band consisting of J"^ = 1~, 3~, • • • due to the 
I30 deformation, which may have a link to the positive-parity bands as a parity doublet. 



C. Energy levels 

In the obtained wave function $amd(-/Vo^''), there exist local minimum solutions in the energy curve as a function 
of A^os- They are also local minima with respect to the deformation /32 and considered to be approximate intrinsic 

wave functions for the corresponding deformed bands. In order to satisfy the orthogonality among the levels and to 
obtain better wave functions, we calculate the level scheme by superposing the spin-parity projected wave functions. 

In the calculations with the case(i) interaction, wc use Nj^t^ =62^+\ 65^+\ 68^+\ 70^+\ 72^+\ 76^+\ 80^+\ 86^+\ 
92^+) with prolate shapes and Not'' = 65^,'^^'' with the oblate shape as the basis for the positive-parity states. For 

the negative-parity states, wc adopt N^^t^ =62'--\ 63(-), 65(-), 66(-), GS^"), 69^^), Tl^"', 65(+), 68(+), 65^^^ The 
obtained level scheme is shown in Fig. 6. By analysing such properties as the dominant components and E2 transition 
strengths of the obtained levels, we classify the levels into groups, which are labeled in Fig. 6(b) and (c). 

In the positive-parity levels, we find K'^ = 0+ rotational bands corresponding to the shape coexistence: (B)the 
oblate band and (A3)the superdeformed band. The states in the groups (Al) and (A2) are built from the mixing 
of the spherical intrinsic state $amd(-^os = 62^"'"^) and the normal prolate state <i>AMD(A^os = 65*^+-'), while the side 
band (A2'(K=2)) of the normal deformation appears from the dominant -Pm2^amd(-^os = 65^+') component due to 
the triaxial nature of the intrinsic state. In Fig. 7, we plot the calculated excitation energies as a function of a spin 
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FIG. 3. The energy curves for the negative-parity states obtained by the constrained AMD calculations. The constraint of 
the total number of harmonic oscillator quanta, (iVos) = A^os, is imposed on the parity projected AMD wave function. The 
energies of the spin-parity projected states, J'mo^'^^amdI-^ob^^)' plotted as a function of Nos in the upper panel(a) and as 
a function of the deformation /32 of ^AMui-^o^^) ^ lower panel(b). The results are obtained with the interaction(i). 
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FIG. 4. Density distributions of the intrinsic states, <5amd(A^3^') for (a) the spherical state with N^t^ = 62'+-', and the 
prolate solutions with N^t^ = (b) 65^+\ (d) 68^+\ (e) 70'+^ (f) 76*+^ (g) 62(-), (h) 65^-', (i) 68^-\ and (c) the oblate 
solution with N^^^ = GS*-^', which are obtained by using the interaction (i). The intrinsic system is projected onto the Y-Z 
plane, and the density is integrated along the X axis, where X, Y and Z axes are chosen as (A^) < (Y'^) < (Z^) and 
(XY) — (YZ) = (ZX) = 0. The deformation parameters p2 are written at the bottom of the figures. The size of the frame 
box is 10 fmxlO fm. 
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FIG. 5. Density distribution of the intrinsic states, $amd(-/Vc;s ) for (a) the the spherical state with 
prolate solutions with iVi^^ = (c) 65^+' and (d) 70(+), (b) the oblate sohition with A'!*^ = 65*+^ and (e) the a-cluster 
state with A^'os'^'' = GS*-"*"^, which are obtained by using the interaction(ii). The intrinsic system is projected onto the Y-Z 
plane, and the density is integrated along the X axis, where X, Y and Z axes are chosen as {X^} < (Y^) < {Z^} and 
(XY) — (YZ) — (ZX) = 0. The deformation parameters /32 are written at the bottom of the figures. The size of the frame 
box is 10 fmxlO fm. 



J(J+ 1). We also show the experimental data of the bands labeled "band 2" and "band 1" in Ref. [12], which are 

assigned to be the normal-deformed band and the superdeformed band, respectively. Since the band (A2) contains the 
configuration mixing, the calculated level structure is not consistent with the rigid rotor levels and does not well agree 
to the experimental data of the normal-deformed band (labeled "band 2" in Rcf. [12]). Moreover, it is higher than 
the superdeformed band in the present results(i), however, we tentatively assign this band as the normal-deformed 
band because of rather strong E2 transitions in (A2) band. On the other hand, the calculated energy spectra of the 
superdeformed band behave as that of the rigid rotor and agrees to the experimental data for "band 1" , though the 
excitation energy of the band-head state is slightly higher than the experimental data. The dominant components 
of the superdeformed band in the low-spin levels is P^Q$Jjy[j-,(iVos = 70^+^), which have a ^^C+^^Si-like clustering. 
Another interesting character of the superdeformed band is its particlc-holc property. As discussed later in detail, 
-Pmo^amd(-^os = 70^+)) has about 90% overlap with the band-head of the superdeformed band, and its intrinsic 
state $amd(.^os = 70^+^) is regarded as a "8p-8/i" state. This particle-hole feature of the superdeformed band is 
supported by the experimental fact that the 0+(5.21 MeV), 2+(5.63 MeV), and 4"'"(6.54 McV) in the superdeformed 
band are strongly populated in 8-nucleon-transfer reactions, ■^^S(^^C,a)"^°Ca [11]. However, it is not consistent with 
the GCM calculations of Skyrme HF-I-BCS [18], where the superdeformed band can not be described by a certain 
np-nh configuration but contains a mixture of various configurations such as Ap-Ah, 6p-6h, and 8p-8h states. In order 
to check the stability of the superdeformed band for the choice of the intrinsic wave functions to be superposed, we 
use another set of the basis as {k} = {N^t^} = {62^+\ 64(+), 65^+), 66^+), 67^+), 68^+), 70^+), 72^+)} and find that 
the level scheme and the E2 transition strengths in the superdeformed band are stable against the choice of the basis, 
and the properties of the states in bands (Al) and (A2) are qualitatively unchanged. In the presetn calculations, the 
number of the basis in the superposition is 10 at most and is much smaller than the GCM calculations in Ref. [18]. 
It also should be noted that, we did not obtain the basis with the quadrupole deformation < /32 < 0.2 as shown in 
Fig. 2. For the detailed investigation of the band mixing, it is expected to be important to perform GCM calculations 
with a large number of the basis concerning the two-dimensional constraints on the A^os and p2- The present results 
of the weak mixings between different np — nh configurations seem to be consistent with the experimental facts, 
however, we comment that the pairing correlations, which were taken into account in Ref. [18], are ignored in the 
present calculations though we think they are partially included through the spin-parity projections and superposition 
of the basis. Since the pairing correlations increase configuration mixings in general, it would be valuable to estimate 
the pairing effects on the stability of the superdeformed band. 

In addition to the prolate bands, the oblate band appears at almost the same excitation energy as that of the 
superdeformed band. The moment of inertia is smaller than the superdeformed band. The oblate deformation is 
predicted also in the macroscopic-microscopic calculations [17] and in the mean-field calculations [18,16]. However, 
since an oblate band has not been experimentally assigned yet, the existence of an oblate shape in '"'Ca is an open 
problem. We also find highly excited bands (A4) and (A5) which arise from the developed ^^C-|-^^Si clustering. The 
details of these bands are discussed in the next section. 

In the negative-parity states(Fig.6c), the lowest 3~ and 5~ states appear from almost the spherical states. The 



9 



theoretical excitation energies are in good agreement with the experimental data of the 3j~(3.74 MeV) and the 
5]^ (4.49 MeV). We find two K"^ — 0" rotational bands (C2) and (C3) which arise from the parity-asymmetric 
shape in the prolate states. Since the energies of the 3^ and 5~ states in these (C2) and (C3) bands arc raised 
by the mixing with the lowest 3~ and 5^ states, the level structure is somehow out of that for the typical K'^ = 
0~ rotational band. The If state(\l/^Z;^) in the band (C2) is dominated by -Pmq^ ^amd(-^os = 65*^~^) with the 
overlap |(*;^=i|-Pmo^^amd(-^os = 65(~)))|^ ~ 90%, while the 1^ state($^~2) in the band (C3) has the dominant 
^mo"'^*^amd(^os = 69'^^^) with 70% overlap. We should note that the If and 1^ states have large overlap with also 
the 1~ states projected from the intrinsic states for the positive-parity, $amd(-^os = 65^+') and $amd(-^os = 68^+'), 
respectively. It is reasonable because $amd(-/Vos~^) is almost the same state as the $amd(-^os^^) when the same 
constraint {Nos) = -^os is imposed. In order to find the possible parity doublet of the superdeformed band, we pay 
special attention to the $amd(-^os = 68^+^) component in the largely deformed negative-parity band (C3) and in the 
superdeformed band (A3) with positive-parity. The overlap |(*I'^=2l-fA^o^^AMD(-^os — 68(+)))p in the bandhead of 
the K'^ = 0~ band(C3) is about 70%, while the bandhead 0+ of the superdeformed band(A3) has 80% overlap with 
-Pmo°^amd(-^os = 68^+'). Then it leads to an interpretation that the largely deformed negative-parity band (C3) and 
the positive-parity superdeformed band (A3) are approximately the parity doublets K'^ = 0~ and K'^ = 0+ arising 
from the intrinsic state <I>amd(A'^os = BS'^^)), which has the parity asymmetric shape with the I30 deformation due to 
the ^^C-|-^*Si-like clustering as shown in Fig. 4. 

Next we discuss the energy levels calculated with the interaction(ii). We use JVi±) =62(+), 65(+), 68(+), 70(+), 72(+) 
with prolate deformations, Not^^ = 65^^^ with the oblate shape, and Not^ = 66^"^^ 68^'*'^ 70i'''^ with a-cluster-like 
structure as the basis in superposing PiiK^AMDi^)- ^he calculated excitation energies are shown in Fig. 8. The 
rotational bands of oblate deformation, prolate one and superdeformation coexist as well as in the results(i). The 
excitation energies of these bands with the intcraction(ii) are much higher than those of the rcsults(i), and also 
are inconsistent with those of the experimentally known low-lying states. For example, the oblate band (B) starts 
from 17 MeV excitation energy, while the bandhead energy of the superdeformed band (A3) is 26 MeV. One of the 
characteristics of the results(ii) is that a a-cluster-like band (D) appears at the relatively lower energy than the 
superdeformed band. 

It has been known that the low-spin states in the normal-deformed band built on the 0^(3.35 MeV) state are strongly 
populated in a-transfer reactions, and they have relatively larger a-spectroscopic factors than other low-lying states 
[4-7]. Therefore, it is expected that the components of the a-clustcr state should be contained in the normal-deformed 
band. However, in the present calculations we failed to dcsc;ribc the a-cluster components in the low-lying states, 
and even in the results(ii), the components of the a-cluster state do not mix with the normal-deformed band (A2) 
except for the 2+ state. It is because the energy difference between the bands (A2) and (D) are large in the present 
calculations. In order to solve this problem, it may be necessary to introduce a suitable nuclear interaction and 
extended model wave functions. 



D. E2 transitions and moments 



We calculate the E2 transition strengths and the electric quadrupole moments for the excited states. We also extract 
quantities related to the intrinsic deformation from the physical observables such as B{E2) and electric quadrupole 
moments. By using the collective model approximation [35], we define the intrinsic quadrupole moment and the 
deformation by the E2 transitions as, 



16n B{E2, A ^Jf) ^^ttQI 
^° Y 5e2 (^i^20|J/i^)2 ' Vl67r3i?2Z' ^ ^ 

where {JiK2Q\J jK) is a Clcbsch-Gordan coefficient, t stands for "transition". We also extract the intrinsic quantities 
from the observable electric quadrupole moments Q as follows. 



l (J + l)(2J + 3) Qs_l^^ ^7) 

Here s stands for "spectroscopic" . 

In the table I, II and III, we show the calculated B{E2) and Q moments, and related quantities. The E2 transitions 
are small among the states in the group (Al), which are dominated by spherical states and considered not to be a 
rotational band of a deformed state. On the other hand, the prolate band (A2) and the oblate one(B) have moderate 
transition strengths, while the E2 transitions are remarkably strong in the superdeformed band (A3). 
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FIG. 6. Excitation energies of the levels up to J = 8 in *°Ca. The experimental data of the levels observed in Ref. [12] are 

shown in (a). Figures (b) and (c) show the theoretical results for the positive-parity and negative-parity states, respectively, 
calculated by the superposition of the AMD wave functions. The interaction(i) is used. The excitation energies of the higher 
spin(J > 8) states will be shown in the next figures. 
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FIG. 7. Excitation energies of the positive-parity states in Ca as a function of a spin J(J + 1). Theoretical values shown 
in the lower panel (b) are those obtained with the interaction (i). The energies are calculated by the superposition of the AMD 
wave functions. The upper panel (a) shows the experimental data for the normal-deformed band and the superdeformed band, 
which are labeled "band 2" and "band 1" in Ref. [12], respectively. 
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FIG. 8. Excitation energies of the positive-parity states in ""^Ca obtained with the interaction(ii). The energies are calculated 
by the superposition of the AMD wave functions. 
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FIG. 9. The deformation parameters /?* extracted from the B{E2) values for the intrarband transitions. The labels of the 
bands correspond to those in Fig.6. The results are for the calculations with the interaction(i). 

In the results with the interaction(i), the deformation parameter /3* extracted from the B{E2) is consistent with /?" 
evaluated with the Q moment except for negative-parity bands and high spin states in the prolate bands (A1),(A2), 
and (A3). In Fig. 9, wc show the behavior of the transition deformation /?* with the increase of the total-angular 
momentum J. In the oblate band (B), the intrinsic deformation fi is almost constant with the increase of J, because 
the oblate states do not mix with the other shapes. On the other hand, /3* in the bands (Al), (A2) and (A3) changes 
as J increases. The deformation of the low-spin states in the superdeformed band is the largest as /3* « 0.5. With the 
increase of J, /3* of the superdeformed band (A3) decreases, while that of the normal-deformed band (A2) enlarges, 
and the /3*-lines for these two bands cross each other around J = 8. We remark that the state mixing occurs in the 
region J > 6, and the components contained in (A2) and (A3) change around J = 8. As a result, the transition from 
the 8+(A2) to 6+ (A3) is significant strong(see Table I). Therefore, another assignment may be possible as S'^{A2), 
10+(A2) and 12+(A2) belong to the superdeformed band. In fact, /?* extracted from the B{E2; J ^ J — 2) among the 
0+{A3), 2+ {A3), 4+ (A3), 6+ (A3), 8+(A2), 10+ (A2) and 12+ (A2) are so large as shown by filled triangles in Fig.9, 
that these states can be assigned to compose a rotational band with the strong E2 transitions. In the experimental 
data by Ideguchi et al. [?], they observed the deviation from the constant moment of inertia in the superdeformed 
band and also the inter-band gamma transitions between normal deformation and supcrdcformation. These facts 
may give an indication of the possible band mixing or some structure change around J ~ 8 in the superdeformation. 
Further fine measurements of E2 transition strengths will be helpful to see the details of the band mixing and to 
establish the assignment for the high spin states of the superdeformed band. 

In the results with the interaction (ii) (table III), we also find strong E2 transitions in the superdeformed band 
(A3). The transitions of the prolate band (A2) are weaker in the results(ii) than in the case(i). As mentioned before, 
the = 0+(D) and K'^ = 2+(D') bands with the a-cluster structure appear below the superdeformed band in the 
results(ii). These bands have strong intra-band E2 transitions, which give the larger transition deformations /3* Ri 0.4 
than the prolate band (A2). 

Let us make a comparison of intra-band transitions between the present results and the experimental data, and 
also compare them with other theoretical calculations. The theoretical and experimental values of Qq for the normal- 
deformed band and the superdeformed band are listed in tabic IV. The present results for the superdeformed band 
are in reasonable agreements with the experimental data. They are also consistent with the results by GCM-I-HFBCS 
calculations. In the shell model calculations by Poves et al., large Qq values were obtained with the truncated 
calculations within 8p-8h configurations by using effective charges = 1.5 and Qi, = 0.5 for protons and neutrons, 
however, the Qq get eroded when other particle-hole configurations are mixed in the calculations. The Qq moments 
of the band (A2) are Qo = 57—112 fm^ and comparable to the experimental values Qq = 74:± 14 fm^ for the normal 
deformation(band 2) in Ref. [12], but are smaller than the experimental data Qq = l.litO.l x 10^ and 1.3±0.2tzmesl0^ 
fm^ in Ref. [36]. Our results are consistent with the GCM-I-IIFBCS calculations except for the 4+ state, which has 
a small Qq in the GCM-I-HFBCS result. As shown in table III, we obtain the a-cluster band (D) with strong E2 
transitions in case of the interaction(ii). The transition moment of the band (D) is Qq « 100 fm^ and well agrees to 
the experimental data in Ref. [36]. Therefore, there might be an alternative assignment that the band (D) corresponds 
to the normal-deformed band built on the 0^(3.35 MeV), though the excitation energy of the band (D) is much higher 
than the experimental data in the present calculations. 

The K = 2 band were experimentally established based on the measurement of B{E2) [1,37]. The observed intra- 
band transitions are B{E2;2,+ 2+) = 730 ± 300 (e^ fm^) and B{E2;A+ 2+) = 200 -t- 150 - 75 (e^ fm^) [1]. 
The properties of this band were studied by Sakuda and Ohkubo with Q!-|-^^Ar(2+)-cluster model [9,10], which well 
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rcpoduccs these intra-band transitions. In the present calculations, we find the K'^ = 2+ band (A2') with the prolate 
deformation in the results (i), and the — 2"^ band (D') with the a-clustcr structure in the case of (ii). In the band 
(A2'), the calculated E2 transition strengths are B{E2:3+ 2+=95 (e^ fm'') and B{E2;A+ 2+=38 (e^ fni''), 
which are smaller than experimental B{E2) values. On the other hands, the E2 transition strengths in the band (D') 
are large enough as B{E2; 3+ 2+) = 390 (e^ fm'') and B{E2; 4+ 2+) = 129 (e^ fm^), which well correspond to 
the data within the experimental error bars. It means that the deformation of the K = 2 band can be described by 
the a-cluster band (D'), which is obtained with the interaction (ii). 

Next we discuss the moments of the negative-parity states. For the negative-parity states in the group (CI), the 
magnetic dipolc moments of the 3~ and 5~ states arc /i(3~) — 1.66 fiN and /i(5~) = 2.76 iin, which well agree to the 
experimental values for the lowest 3]~(3.74 MeV) and the lowest 5J~(4.49 MeV); A*(3^) = 1.65 /xjv and At(5~) = 2.6 /Ltjv- 
In the prolate bands (C2) and (C3), the E2 transitions are fragmented among the bands because the 3~ and 5~ states 
of these bands contain configuration mixing due to the existence of the lowest 3~ and 5~ states in the group (CI) 
(see table II). Then, we estimate the intrinsic deformation by looking into extracted from the Q moments instead 
of (3*. According to the calculated /S" of the band-head state(l^), the deformation of the band (C3) is evaluated as 
= 0.39. This is smaller than that of the superdeformed band (A3), which we raughly assign to the parity doublet 
of (C3). 

IV. DISCUSSION 
A. Molecular bands 

Appearence of molecular resonances is one of the important cluster aspects in sd-shell nuclei. In the large prolate 
deformations of ^°Ca, we find the trend of developed ^^C+^^Si clustering which leads to the superdeformed band as 
mentioned in the previous section. Therefore, it is natural to expect possible existence of ^^C+^^Si molecular bands 
and their link with the superdeformed band in ^°Ca. 

It has been known that molecular resonances appear in such systems as ^^C-t-^^C, ^^C+^^0 or ^^0-|-^^0 [38,39]. In 
the microscopic calculations, there have been theoretical efforts to connect the molecular resonances with the low-lying 
deformed bands [22,40-42]. We should stress that the recent AMD studies [22,42] are the first theoretical works which 
can microscopically describe both the coexisting deformed states in low-energy region and the molecular resonances 
in the high-energy region without relying on assumptions of constituent clusters. In Ref. [22], the relation between 
the theoretically predicted superdeformed band in ^^S and the ^^0+^®0 molecular bands has been studied. It has 
been suggested that the superdeformed band with a considerable amount of the ^^O+^^O component is regarded as 
the lowest nodal band in a series of molecular bands, while the excited ^^0-|-^^0 bands arise due to the 

inter-cluster excitation and correspond to the observed molecular resonances. 

In the present results of the case(i), we find the excited rotational bands (A4) and (A5) as shown in Fig. 7. The 
former and the latter exist at about 15 MeV and 25 MeV higher excitation energies than that of the superdeformed 
band, respectively. These bands appear due to the superposition of the large prolate deformations with the ^^C-|-^^Si 
cluster structure. In Fig.lO, we show the overlap of the states(^';[+) in the prolate bands (Al), (A2), (A3), (A4) and 
(A5) with a single AMD state: |(^';{+|^'mo*amd(^o^'')) The low-spin(J < 4) states in the superdeformed band 
is dominated by ^Mo*i'AMD(-'^os ~ 70^+') with about 90% overlap. We should stress again that the ^^C+^^Si-like 
clustering actually appears in the <i>AMD(A^os — 70^+^) (see Fig. 4). On the other hand, the states in the excited band 
(A4) contain the AMD states with larger A'os, and the components in the band (A5) shift toward the further large Nos 
region. Since the di-nuclear clustering develops further and two clusters go away from each other with the increase 
of Nos, these excited bands (A4) and (A5) should be described by the excitation of inter-cluster motion. In fact, we 
calculate the oscillator quanta of the subsystem consisting of the nucleons within each cluster(^^C or ^^Si), and find 
that each of the cluster is almost written by OHlo configurations. In other words, the increase of A^os corresponds to 
just enlarging the inter-cluster distance, and hence, the rotational bands (A4) and (A5) arise due to the excitation 
of relative motion between ^^C and ^*^Si clusters. Moreover, it is suggested that the bands (A4) and (A5) may have 
tails in the inter-cluster motion because the components of these bands are spread into the broader A'os ranges than 
that of the superdeformed band. These results suggest an interpretation that the superdeformed band (A3) and the 
higher bands (A4) and (A5) are regarded as a series of '^^C-l-^^Si molecular bands. Namely, the bands (A4) and (A5) 
appear as the higher nodal states built above the the superdeformed band. These results have a good analogy to the 
feature of the superdeformed band and the -^^0-1-^^0 bands in ^^S discussed in Ref. [22]. 

In the experimental side, the ^^C-|-^*Si resonances have been observed in the elastic scattering data in the backward- 
angle region [24,38]. The resonances at the energy Ecm = 26.0 and 30.2 MeV, which correspond to the excited 
states of '^''Ca at £J = 39 and 44 MeV, are assigned to be J = 18 states from the angular distributions. There 
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TABLE I. Electric quadrupole moments Q (e fm^) and B{E2) (e^fm'*) for the positive-parity states calculated with the 
interaction(i). The extracted intrinsic quadrupole moments (fm^) and deformation parameters are also listed. The definitions 
of Qq" and (3*'" are explained in the text. The labels of the bands correspond those in Fig.6. 



tarisition to 



J" 


band 


Q 


Qo 


13' 


final 


band 


B{E2) 


Qo 




positive-parity bands 


2+ 


(Al) 


-10 


34 


0.13 


0+ 


(A2) 


12 


24 


0.10 


4+ 


(Al) 


-10 


28 


0.11 


2+ 


(Al) 


28 


31 


0.12 


6+ 


(Al) 


-13 


33 


0.13 


4+ 


(Al) 


27 


29 


0.11 


8+ 


(Al) 


—15 


36 


0.14 


6^ 


(Al) 


38 


34 


0.13 


10+ 


(Al) 


-22 


50 


0.19 


8+ 


(Al) 


54 


40 


0.16 


12+ 


(Al) 


-30 


68 


0.27 


10+ 


(Al) 


96 


53 


0.21 


prolate 
























2+ 


(A2) 


-17 


61 


0.24 


0+ 


(A2) 


64 


57 


0.22 


4+ 


(A2) 


-24 


65 


0.26 


2+ 


(A2) 


116 


64 


0.25 


6+ 


(A2) 


-29 


72 


0.28 


4+ 


(A2) 


128 


64 


0.25 


8+ 


(A2) 


—43 


102 


0.40 


6^ 


(A2) 


200 


78 


0.31 


8+ 


(A2) 








6+ 


(A3) 


161 


70 


0.27 


10+ 


(A2) 


-50 


116 


0.45 


0-4- 

8^ 


(A2) 


370 


105 


0.41 


12+ 


(A2) 


-56 


126 


0.49 


10^ 


(A2) 


433 


112 


0.44 


superdeformation 




















2+ 


(A3) 


-37 


131 


0.51 


0+ 


(A3) 


330 


129 


0.51 


4+ 


(A3) 


-47 


129 


0.51 


2+ 


(A3) 


476 


129 


0.51 


6+ 


(A3) 


-47 


118 


0.46 


A 4- 

4+ 


(A3) 


455 


121 


0.47 


8+ 


(A3) 


-40 


95 


0.37 


6+ 


(A3) 


309 


97 


0.38 


10+ 


(A3) 


-36 


83 


0.33 


8+ 


(A3) 


222 


81 


0.32 


12+ 


(A3) 


-43 


97 


0.38 


10+ 


(A3) 


210 


78 


0.31 


oblate 
























2+ 


(B) 


19 


—67 


—0.26 


0^ 


(B) 


91 


—68 


—0.26 


4+ 


(B) 


24 


-66 


-0.26 


2+ 


(B) 


128 


-67 


-0.26 


6+ 


(B) 


25 


-63 


-0.25 


4+ 


(B) 


139 


-67 


-0.26 


8+ 


(B) 


25 


-60 


-0.24 


6+ 


(B) 


143 


-66 


-0.26 


in+ 


(B) 


25 


-57 


-n.22 


8+ 


(B) 


144 


-65 


-0.26 


12' 


(B) 


2:i 


-52 


-0.21 


10 


(B) 


iii 


-()5 


-0.25 


K = 2h 


and 






















2+ 


(A2' 


K=2) 


13 


47 


0.18 














3+ 


(A2' 


K=2) 








2+ 


(A2' 


K=2) 


95 


52 


0.20 


4+ 


(A2' 


K=2) 


-9 


61 


0.24 


2+ 


(A2' 


K=2) 


38 


57 


0.22 


5+ 


(A2' 


K=2) 


-12 


53 


0.21 


3+ 


(A2' 


K=2) 


25 


36 


0.14 


6+ 


(A2' 


K=2) 


-18 


62 


0.24 


4+ 


(A2' 


K=2) 


83 


60 


0.23 


7+ 


(A2' 


K=2) 


-20 


63 


0.25 


5+ 


(A2' 


K=2) 


58 


47 


0.18 


8+ 


(A2' 


K=2) 


-22 


63 


0.25 


6+ 


(A2' 


K=2) 


104 


61 


0.24 


9+ 


(A2' 


K=2) 


-24 


64 


0.25 


7+ 


(A2' 


K=2) 


113 


62 


0.24 


10+ 


(A2' 


K=2) 


-22 


57 


0.22 


8+ 


(A2' 


K=2) 


96 


56 


0.22 


11+ 


(A2' 


K=2) 


-25 


63 


0.25 


9+ 


(A2' 


K=2) 


114 


60 


0.24 


12+ 


(A2' 


K=2) 


-28 


67 


0.26 


10+ 


(A2' 


K=2) 


99 


55 


0.22 
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TABLE II. Electric quadrupole moments Q (e fm^) and B{E2) (e^fm*) for the negative-parity states calculated with the 
interaction(i). The extracted intrinsic quadrupole moments (fm^) and deformation parameters are also listed. The definitions 
of Qq'' and (3*'^ are explained in the text. The labels for the bands are the same as those in Fig.6. 
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(C2) 


150 
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are other experimental implications of ^^C-|-^^Si molecular states in the low-energy fusion cross section [23], where 
distinct structures have been found in the energy region _Ecm=20 — 30 MeV above the ^^C-|-^*Si threshold. These 
are the candidates of the molecular resonances which might correspond to the -'^^C+^^Si molecular bands (A4) and 
(A5) obtained in the present results. In order to make further investigations of the molecular bands and give more 
quantitative discussions, it should be important to superpose a large number of the AMD states. It also may be 
effective to take into account the excitation inside the clusters as coupled channel calculations for the description of 
the detailed band structure. 



B. single- particle orbits 

In an AMD wave function, the single-particle wave packets, ipi, in Eq.(3) are not orthogonal to each other. In order 
to study the mean-field character, it is useful to transform them to the HF-like single-particle orbits {'p^^}^ which 
are orthonormal to each other and form the total wave function equivalent to the original AMD wave function. We 
extracted the single-particle orbits {(/'a^} from <&amd(A"os^''), as explained in Refs [32,43], and analyze them. 

As mentioned before, the low-spin states of the superdeformed band is dominated by the AMD state, 
-Pmo^amd(-^os = TO^"*"^). We find that the intrinsic state, $amd(A^os = 70^+^), contains eight /p-like single-particle 
levels, which are the highest 4 levels among proton orbits and the highest 4 levels among the neutron orbits. Each of 
these levels has dominant parity-odd component as the ratio P{—) ^ 95%. The density distributions of these levels 
are shown in Fig. 11(a) and (b). Hereafter, we use the label a(=l,- • -,20) of {f^^} for the oth highest single-particle 
levels in the proton orbits. It is found that the properties of the eight valence levels in $amd(-/Vos = 70^+^) are 
associated with the longitudinal fp orbits which have l3,±i and ¥3,0 angular dependent terms. Therefore, we regard 
the superdeformed state as a "Sp-Sft," state, though it is not equivalent to pure particle-hole states in the mean-field 
picture. The origin of the parity asymmetric shape with ^^C-|-^^Si cluster in i'AMD(Afos = 70^+^) will be discussed 
later. 

The prolate band (A2) contains the significant component of the AMD state projected from the intrinsic state, 
^amd(-^os = 65^+'). The valence orbits in this state are the eight levels which contain 50% parity-odd components 
and 50% parity-even ones as shown in figures 11(d) and (e). It means that the parity projected state ^amd(-^os = 
65*^+)) can be approximately described by the mixing of 2p-2h, ip-Ah and 6p-6h configurations. We comment that 
^amd(A^os = 65*^+-') is not a local minimum in the energy curve as a function of iVos nor /3. 

The oblate solution $amd(A^os = 65^^^) does not mix with the prolate states, and independently forms the oblate 
band (B). This state can be convincingly described by a Ap-Ah state because 4 valence nucleons occupy almost parity- 



16 



TABLE III. Electric quadrupolc moments Q (e fm'^) and B{E2) (e'^fm*) calculated with the interaction(ii). The extracted 
intrinsic quadrupole moments (fm^) and deformation parameters are also listed. The definitions of Qo " and p*'" are explained 
in the text. The labels for the bands correspond to those in Fig.8. 
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oscillator quanta 

FIG. 10. The overlap between the calculated excited states '^i,'^ and the prolate AMD states Pj^o*amd(-^os^')- values 
IPiio^AMDi^ot^ ))\'^ fo'" states in the bands (Al), (A2), (A3), (A4) or (A5) are shown by open circles, open triangles, 
open squares, filled circles or filled triangles, respectively. 
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TABLE IV. Intrinsic quadrupole moments (fm'^) of the normal-deformed band and the superdeformed band extracted from 
B{E2) with the definition in Eq.(16). The present AMD results for the normal-deformed band and the superdeformed band 
are those for the bands (A2) and (A3) with the interaction(i), respectively. The values in parensis are those for the assignment 
of the 8^, 10^, 12+ states in the band (A2) to the superdeformed states. The shell model calculations by Poves et al. [19], and 
the GCM-I-HFBCS calculations by Bender et al. [18] are also listed. The experimental data are those taken from B{E2) data 
in Ref [36], and Qq values determined by Doppler-shift analysis in Refs. [12] and [13]. "For the values of Refs. [12] and [13], a 
single Qo value is assumed for the entire transitions in each of the "band 2" or the "band 1" . 
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odd orbits with angular parts related to Vs.is- The property of a valence single-particle level of the four is shown in 

Fig. 11(c). The particle- hole nature is foimd also in the a-cluster states $amd(./Vos = 66^"^^) and ^amt>{Nos = 68^^^) 
obtained by the interaction(ii). Namely, these states have a Ap-Ah feature as four valence neutrons occupy /p-like 

orbits with 75% negative-parity component in the $amd(-^os = 68^^'') as seen in Fig. 11(f). 

Thus, the single-particle levels in the intrinsic states disclose the many- particle many- hole characters in ''"Ca. On 
the other hand, the cluster aspect has been found in prolate states which have parity-asymmetric shapes. It raises a 
question how the cluster aspect links with the mean- field aspect. Here we examine the cluster aspect while paying 
attention to parity asymmetry of the single-particle orbits. In Fig. 12, we show the ratio of negative-parity component 

of the single-particle proton orbits in the prolate states <I>amd {Not^ ) . The behavior of the neutron orbits is almost the 
same as that of the proton orbits. At almost the spherical Nog = 62 region, 12 protons(i^^Z;^_;^2) occupy almost parity- 
even orbits and 6 protons (</?^£i3_ig) are in almost parity-odd orbits, which indicates the approximately srf-shell closed 
configuration of this state. With the increase of A'os, the parity-odd components in the highest four leveh((p^^-^^_^) 
increase and these four levels become almost pure parity-odd orbits at A^os = 70. It means that four protons in the 
srf-shell are excited into the /p-shell. In the state ^AMu{Not^ = 70^+^), each of the higher ten levels ip^Ei_iQ has 
approximately good parity. Instead, the parity symmetry breaking in single-particle orbits is rather remarkable in 
the 11th and 12th \evels{(p^Eii which are originally the lowest sd-orbits as shown in Fig. 12 (b). In order to see 

the effect of the single-particle features on the collective shape of the total system ^AMT){Not^ = 70^+^), we show the 
density distribution summed up for the protons in the lower \evels{ip^Eii_2o) and that for the protons in the higher 
levels(</3Q£i_iQ) in Fig. 13. The most striking thing is that 10 protons in the lower levels form a parity asymmetric 
shape, which reminds us of the clustering in ^°Ne system. These results imply that the parity asymmetry in 

this state dominantly originates not in the valence nuclcons near the Fermi surface but in the a-cluster-likc correlation 
in the core part. It is inconsistent with an mean- field picture that nucleons at the surface may contribute such an exotic 

shape as the octupole deformation. We here remind the reader that this state($AMD(.^os^^ = 70^+')) is the dominant 
component of the superdeformation and has the ^^C+^^Si-cluster structure(Fig. 13(c)). The ^^C+^^Si clustering is 
also seen in the aspect of molecular resonances. This seems to be inconsistent with possible ^^S-|-2a clustering, which 
may be naively expected from the 8p-8h nature of the superdeformed state. By the analysis of single-particle orbits, 
we can interpret the appearance of the ^^C cluster as follows. The origin of the -"^^C cluster in the superdeformation 
is understood by the correlation within 12 nucleons; 4 sd-shell nucleons with the a correlation and the 8 valence 
nucleons in the /p-like orbits. It leads to a new picture of the coexisting cluster and mean- field aspect. Namely, 
the deeply bound nucleons compose the parity-asymmetric core ^°Ne with the a-cluster correlation, while the single- 
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FIG. 11. (Color online) jDensity distributions of the HF-like single-particle orbits (p™ for protons in $AMD(Afos )• The 
results are for the calculations with the interaction(i) except for the panel(f). The first(iy3™i) and the third(<^J^£3) highest proton 
orbits in the the dominant component("I>AMD(Afos'^' ~ 70*^^') of the superdeformation are shown in (a) and (b), respectively, 
(c) The highest orbit, ip^Ki, in the oblate solution with Not'' = 65'"''''. (d)The first and (e)the third(</3a£3) orbits in the 
prolate solutions with N^t^ = Qb'^^K (f) The highest orbit in the Q-cluster state with A'^s'^' — GS'"*"' obtained in the 
results with the interaction(ii). In the left side, the intrinsic state is projected onto the Y-Z plane, and the density is integrated 
along the X axis, where X, Y and Z axes are chosen as (X^) < (F^) < (Z^) and (XY) = (YZ) = (ZX) = 0. The right figures 
are for the surface cut of the density. The percentage P{ — ) of the negative-parity component in each single-particle orbit is 
also listed. The box size is 10 fm. 

particle levels at the energy surface have particle-hole nature. In the further large A'os region, the particlc-holc nature 
disappears because of the spatial development of ^^C+^^Si clustering. 

V. SUMMARY 

We studied deformed states in "^^Ca with antisymmetrized molecular dynamics(AMD) by using effective nuclear 
interactions which contain finite-range 2-body and 3-body forces. In the framework of a constraint AMD method, 
we performed parity projection before energy variation, and total-angular-momcntum projection after the energy 
variation. The obtained AMD wave functions are superposed to get better wave functions. 

By analyzing the calculated results of the intrinsic states, level scheme, and E2 transition strengths, it was found 
that the rotational bands with various deformations appear in ^°Ca as well as the spherical ground state. Namely, 
above the spherical ground state, the rotational bands arise from the normal-deformed state and the superdeformed 
state as well as the oblate state. In the results with a weaker spin-orbit force, we also obtained the a-cluster-like 
bands. The present results of the superdeformed band reasonably reproduce the properties of the experimental 
superdeformed band build on the Og (5.21 MeV). On the other hand, we could not obtain a satisfactory description 
of the experimental deformed band build on the 0^(3.35 MeV). Possible assignment of this band is to the theoretical 
normal-deformed band, the a-cluster band, or admixture of them. 
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FIG. 12. The percentage P( — ) of the negative-parity component in each HF-hke single-particle orbit for the protons in the 
prolate states $AMD(A/i^^) obtained with the interaction(i). The values P(— ) for the highest ten levels(a = 1 — 10) are shown 
in the upper panel, and those for the orbits(o = 11 — 18) from the 11th to the 18th level are shown in the lower panel. 




(b) 



Nos=70(+) 
a= 11-20 





FIG. 13. Density distributions summed up (a) for the protons in the higher single-paxticle levels {'fia=i-io) and (b) for the 
protons in the lower single-particle levels 20) in the intrinsic state of the superdeformation:$AMD(iVoB = 70^"*"^). The 
total density is also shown in (c). The densities in (a) and (b) are normalized by a factor 4. 
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In the anslysis of single-particle orbits, we found the particle-hole nature of the superdeformed state, the oblate 
state and the a-cluster state, which are dominated by the 8p-8h, Ap-Ah and Ap-Ah configurations, respectively. On the 
other hand, the normal deformation in the present results contains mixing of 2p-2h, Ap-Ah and 6p-6h configurations. 

One of new findings in the present study is that the superdeformed state has a parity asymmetric shape with 
i2(-<_|_28gj_jj]jg cluster structure. From the point of view of a single-particle picture, we discussed the relation between 
the cluster aspect and the many-particle many- hole aspect. We found that the origin of the ^^C clustering in the 
superdeformed state is understood by the correlation of 12 nucleons; 4 sd-shell nucleons with the a-cluster correlation 
and the 8 valence nucleons in the /p-like orbits. It reveals the coexistence of cluster and mean-field aspect in 
the superdeformed state. Namely, the a-clustcr correlation in deeply bound nucleons composing the core ^"^Nc 
plays an important role in the parity asymmetric shape, while the particle-hole nature arises in the single-particle 
levels at the energy surface. The present result of the dominant 8p-8h feature in the superdeformation seems to 
be consistent with the strong population of the members in the superdeformed band observed in 8-nucleon-transfcr 
reactions;'^^S(^^C,Q!)"'"Ca [11]. However, it is inconsistent with the GCM-I-HFBCS calculations [18], where neither the 
Ap-Ah nor the 8p-8h HF solutions can be assigned to the bandheads of the normal-deformed band or superdeformed 
band. One of the advantages of the present framework is the parity projection before energy variation, which is 
considered to be important in arising of the parity-asymmetric deformation. We think that one of the key properties 
which stabilize the "8p-8/i" character in the superdeformed band is the ^^C+^^Si clustering which causes a parity 
asymmetric shape. Compared with the GCM-I-HFBCS calculations [18], the higher correlations are ignored in the 
present calculations. For example, the number of the basis in the superposition is ten at most, which is much smaller 
than the case of the GCM+HFBCS calculations [18]. Moreover, since we adopt the constraint on the total number of 
HO quanta in the energy variation, some basis are missed in terms of GCM with respect to the generator coordinate 
/32, which are often used in the GCM-I-HF calculations. The pairing correlations are also ignored. It is an remaining 
problem to check how the properties of the superdeformation are affected by introducing these effects. 

Above the superdeformed band, we found possible higher rotational bands based on ^^C+^'^Si-clustering in the 
large prolate deformations. It is predicted that ^^C+^*Si molecular bands may be built due to the excitation of the 
inter-cluster motion. The present results suggest that the superdeformed band and these higher molecular bands are 
regarded as a series of ^^C-|-^*Si molecular bands. This means that the ^^C-|-^^Si cluster aspect is rather prominent 
in the superdeformed state than •^^S-|-2a aspect, although the latter may be naively expected from the 8p-8h feature 
of this state. We also suggest the negative-parity bands caused by the parity asymmetric deformation. 

In the present work, we use the stronger spin-orbit force than the nuclear interaction adopted in Ref. [30] to 
quantitatively reproduce the excitation energies of the superdeformed band and negative parity states. With the 
original weak spin-orbit force, the excitation energies of all the excited states in '^''Ca are overestimated. We would 
like to stress again that the intra-band properties such as energy spectra, E2 transitions and intrinsic structure of 
the superdeformed band are not sensitive to choice of either interaction. In order to solve the remaining problems, it 
may be necessary to introduce a suitable nuclear interaction and extended model wave functions. 
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